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1. Introduction 

Since STAPHYS 21 our understanding of the statistical properties of the one-dimensional 
KPZ equation (Kadar, Parisi, and Zhang [1]) for surface growth has improved substantially. 
Thus an obvious scheme would be to sketch the main results and to list the relevant articles. 
While the latter will be done anyhow, it seems to me that a more interesting strategy is to 
explain one central method through which these advances became possible, namely the 
mapping to line ensembles. Given the appropriate line ensemble, to obtain properties of 
physical interest still requires a tough asymptotic analysis for which I refer to the articles 
[2-10] . The link to multi-matrix models is illustrated in our contribution to the Proceedings 
of thelCMP 14 [11]. 

It seems difficult to handle directly the continuum KPZ equation and one has to rely on 
suitable discretizations. In one dimension the slope is locally conserved. Thus a natural 
discretization is to assume that the slope takes only two values and that the dynamics 
is defind through a suitable stochastic exchange rule. Thereby one arrives at the totally 
asymmetric simple exclusion process (TASEP), which is also a standard model for driven 
lattice gases in one dimension [12]. We briefly recall the definition. At each site j G Z 
there is at most one particle. Thus the occupation variables rjj (t) at time t take only the 
values (site j empty) and 1 (site j occupied). Under the stochastic dynamics a particle 
with an empty right neighbor site jumps independently to that site after an exponentially 
distributed waiting time with rate 1 . The master equation for the Markov chain thus reads 



ift(v) = Lf t (v), t>0, 



(1) 
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with generator 

Lf(v) = E%(! - Vj + i)(f(rf' j+1 ) f{v)) (2) 

as acting on functions / over configuration space. Here rf'^ 1 denotes the configuration r/ 
with the occupations at sites j and j + 1 interchanged. 

In the interpretation of surface growth, the height function hj (t) at time t has the height 
differences hj+i(t) — hj(t) = 1 — 2rjj + i(t) with the convention that ho(t) = 2N t , N t 
being the number of jumps through the bond from to 1 in the time span [0,i]. Then 
our dynamical rule simply means that local minima of hj (t) are filled independently af- 
ter an exponentially distributed waiting time. The stationary measures of the TASEP are 
Bernoulli, Pmb p (i]j = 1) = p. Therefore {r]ir]j)p — p 2 = — p) and the stationary 

current is j(p) = p(l — p). Since j"{p) ^ 0, the TASEP is in the KPZ universality class 
for every density p, p ^ 0, 1. More details can be found in [13]. 

For the initial conditions there are three prototypical choices. In fact, while the dynam- 
ical scaling exponent is always z = 3/2, the scaling functions do depend on the initial 
conditions. 

-flat initial data, which is the standard choice for numerical simulations. For the TASEP 
it means 77^ (0) = (1 + (— l) J ')/2. Then (ho(t)) = \t. From a statistical mechanics point 
of view the object of interest are random fluctuations. Their behavior is only partially 
understood. Assuming universality, there is strong evidence for 

ha(t) = ^-£goe* 1/3 fort ^00 (3) 

with the random amplitude £goe distributed as the largest eigenvalue of a N x N GOE 
random matrix in the limit N — > 00 [14,6]. This distribution is known as Tracy-Widom 
at (3 = 1 [15]. Note that instead of centering we subtracted the asymptotic mean. In 
particular, it turns out that (£goe) < 0. As a stronger statistical property, under diffusive 
scaling with scale parameter t 2 / 3 and in the limit t — > 00, the process 

x h-> t~ 1 / 3 (h[ t 2/3 x ^(t) — ^tj , [-J integer part , (4) 

should have the same fluctuations as the top line in Dyson's Brownian motion at (3 = 1 in 
the limit — > 00, see [6] for a complete discussion. 

- droplet geometry, which takes rjj(0) = 1 for j < and r)j(O) = for j > 0. In terms 
of the height, hj(0) = \j\. The asymptotic shape of the droplet is h^ xt j (t) = t(x 2 + l)/2, 
\x\ < 1, which has nonzero curvature in contrast to the flat initial conditions. The droplet 
geometry occurs naturally in the growth of an Eden cluster and other growth processes in 
the plane starting from a point seed. In our context the droplet geometry derives its fame 
from the work of K. Johansson [2], who proved for the height at the origin 

ho(t) = ^-Ccue* 173 for r-> 00 (5) 

with the random amplitude £gue distributed as the largest eigenvalue of a N x N GUE 
random matrix in the limit N — > 00. For the droplet geometry even stronger property (4) is 
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available. The limit is the stationary Airy process [16,3]. (For the sake of notational clarity 
we omit numerical multiplicative factors. They are provided in the references given.) 

- stationary initial conditions. This means that at time t = the r/j (0) are independent with 
Prob p (?7j(0) = 1) = p. The stochastic process i]j(t) is stationary in both variables and 
can be viewed as a two-dimensional field theory in infinite volume. The natural statistical 
object is then the two-point function 

S(j,t) = ( Vj (*)%(0)> P -P 2 . (6) 
It has the scaling form, introducing the compressibility x(p) = p(l — p)> 

S(j,t) - x(4x 1/3 t 2/3 )- 1 ^"((4x 1/3 t 2/3 )- 1 (.7 - (1 - 2p)t)) (7) 

for j — (1 - 2p)t = 0(t 2 / 3 ) and t — > oo. A plot of the scaling function g" can be found 
in [17]. For the TASEP a proof based on line ensembles is under construction [18]. 



2. Nonintersecting Line Ensembles 

Line ensembles are familiar from bosonic and fermionic paths integrals. In our context the 
relevant example are free fermions on the one-dimensional lattice with Hamiltonian 

H = («*0>0" + 1) + + • (8) 

3 

Here the Fermi field a(j) satisfies the anticommutation relations {a(i),a(j)} = 0, 
{a(i),a(j)*} — Sij. Let us denote by xj(s), j = 1, . . . , N, a collection of independent 
symmetric nearest neighbor random walks on Z with jump rate 1 . We condition them not 
to intersect, indicated by NC for non-crossing constraint. Then the A^-particle transition 
amplitude is expressed as 

(xi, . . . ,x N \e~ sH \yi, ... ,y N ) = 

Prob NC (xi(s) =xi,... ,x N (s) = x N \x 1 (0) = y u ... ,x N (0) = y N ) , (9) 

i.e. the random walks start at y\ < . . . < yjv at time and end up at xi < ... < xn 
at time s constrained not to intersect. Since the constraint is purely geometric, it is also 
called entropic repulsion. The transition probability of the TASEP, e Lt , cannot be written 
as a line ensemble with purely entropic interaction. Nevertheless, in a much more hidden 
way there is such a line ensemble, but it depends in a delicate way on the initial conditions 
as will be explained in the following section. 

There is a second example which is formulated rather differently, but close to (8) in fact. 
Following Dyson [19], we consider the stochastic process s \— ► A(s) which takes values in 
the N x N Hermitian matrices and is governed by the linear Langevin equation 

±A(s) = ~A(s) + W(s). (10) 

Here W(s) is an N x N Hermitian matrix whose matrix elements are independent Brow- 
nian motions. We assume A{s) to be stationary in time. A(s) has the real eigenvalues 
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A_jv + i(s) < . . . < Ao(s). They never cross. At a fixed time s, A(s) has the Gaussian 
distribution Z~ x exp[— jjtrA 2 ]dA, known as GUE of random matrices [20]. The eigen- 
values Xj (s), j = —N + 1, . . . ,0, form a line ensemble with purely entropic repulsion. In 
particular, the 2-dimensional random field 

o 

<t> N (x,a) = 2 6(\ 3 (s)-x) (11) 

j=-JV+l 

has determinantal moments constructed out of the extended Hermite kernel [21]. 

As to be detailed, very roughly, the height function in the growth process has the same 
statistics as the top eigenvalue in Dyson's Brownian motion (10), where space x corre- 
sponds to the Brownian motion time s and the growth time t corresponds to the matrix 
size N. This correspondance explains why the universal distribution functions for the KPZ 
equation are expressed through quantities appearing also in the edge scaling of random 
matrices. 



3. Line Ensemble for the TASEP 



We pick the case of step initial conditions, which is the easiest one to explain and well 
illustrates the principle. Thus initially rjj = 1 for j < and rjj = for j > 0. Let us 
label the particle initially at site j by the index — j + 1. We define w(i,j) as the i-th jump 
time of particle number j. Thus w(i,j), i, j = 1, 2, . . . , is an infinite array of independent 
exponentially distributed random variables, Prob(w(i, j) 6 [y, y + dy]) = e~ v dy, y > 0. 
Given the w(i,j)'s one constructs from them a growth process by filling the square 
after the random time w(i,j) starting from the moment when the two adjacent squares 
(i — and (i, j — 1) are already filled. Let us denote by G(m, n) the time when the 
square (m, n) is first filled. Then the cluster grown at time t, {(m, n)\G(m, n) < t}, has 
as its border the height function hj(t) up to a rotation by ir/4. 

G(m, n) can be considered as a sort of dual growth process. It satisfies the random 
recursion relation 

G(m, n) = max{G(m — 1, n),G(m, n — 1)} + w(m, n) , (12) 

G(0, n) = = G(m, 0). G(m, n) has also the interpretation of the optimal time for a last 
passage directed percolation, which can be seen by introducing the directed polymer, lo. It 
is an up/right path on (Z + ) 2 starting at (1,1). The "length" of a polymer is given by 

r(w)= J2 w ^^- ( 13 > 

G(m, n) is the length of the optimal path at fixed endpoints, 

G(m,n) — max t(u>) . (14) 

Since hj(t) is increasing linearly in t, so is G(n + k, n — k) in n, and the scaling behavior 
of hj(t) for large t is determined by the scaling behavior of G(n + k,n — k) for large n. 
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Figure 1. Line ensemble for the TASEP with r = 3. The lines visualize 
the non-crossing constraint. 

For example, G(n, n) is the time it takes for ho(t) to reach the level 2n. Thus the assertion 
(5) follows from a corresponding one for G(n, n) in the limit of large n. 

We are now in the position to construct out of G(m, n) an ensemble of lines with purely 
entropic repulsion. We fix r and introduce the lines hg(j), £ = 0, . . . , — r + 1. hf (j) E R+ 
and he(2j) - h e (2j - 1) > while h e (2j + 1) - h e (2j) < 0. Furthermore h e (j) = for 
\j\ > 2t — 21 — 1, see Figure 1. Each jump of absolute size <5 carries the weight e _<5 / 2 , 
5 > 0. The weight of the whole line ensemble is then the product of the weights of all 
jumps. Finally one has to impose the non-crossing constraint. It is most easily formulated 
through extending j hg(j) to piecewise constant lines in the plane as 

ht(s) = h e (j) for j - 1 < s < j + 1 . (15) 
Note that for each constant piece both endpoints are included. We then require 

(NC) ht-^a) < h e (a) (16) 

for all £, all s, and all endpoints, compare with Figure 1 . 

The idea of the line ensemble is that the top line, ho(j), records the statistics of G(m, n), 
while the others are introduced only for book-keeping purposes. This can be seen by 
generating the line ensemble through the Robinson-Schensted-Knuth algorithm [22] (it is 
the Knuth of TeX). r is regarded as growth time parameter and we set w(l, 1) = h (0, t = 
1), see Figure 2. The up-step (at s = — \) is shifted one unit to the left and the down-step 

(at s = 5) is shifted one unit to the right, yielding ho(j, t = 1.5), j = —1, 0, 1. Then we 
add mass according to h (— 1, r = 1.5) = h (—l,r = 1.5) + w(l, 2), h (0, r = 1.5) = 
ho(0, t = 1.5), h (l, t — 1.5) = ho(l, t = 1.5) + w(2, 1), etc.. In this way the top line 
h (s, t), s G K, r = 1, 2, . . . , satisfies with the identities 

G(t + j, r-j) = h (2j, r) , G(r + j, r - 3 - 1) = h (2j + 1, r) . (17) 

The line h-\ is generated from ho through a fully deterministic rule: When moving the 
up-steps and down-steps of ho to ho there will be excess mass at coalescing steps. This 




x=l x=1.5 x=2 

Figure 2. The RSK algorithm for the TASEP. Light gret areas are mass de- 
posited according to the growth rule. The shaded area is the mass annihila- 
ted in line and copied to line — 1. The lines visualize the updating rule. 



excess mass is recorded as nucleation for h-\. The up-steps and down-steps of h-i move 
deterministically in the same fashion as explained for ho- The nucleation for the line h-2 
is generated by the excess mass of h-\, etc.. Note that the only random input is to the top 
line ho. Of course, it requires an argument, see [10,18], that at growth time r the lines have 
exactly the probability distribution as claimed, namely weighing the jumps exponentially 
and imposing the non-crossing constraint. 

We have to issue a word of caution. For the multi-line growth process probabilities 
referring to two different growth times are meaningful. In particular we may want to 
enquire on ho(s, r) for two different r's. This information is no retraceable from the static, 
fixed r, line ensemble. In case of the TASEP, line ensembles can be used to analyze the 
statistics of hj(t), j G Z, jointly with prescribed initial conditions but not for, say, ho(ti) 
and hoifz) jointly. 

As in (1 1) for Dyson's Brownian motion the line ensemble ht{j, t) determines the ran- 
dom field (p T (x,j) through 

o 

Mx,j)= E 8{h t (j,T)-x). (18) 
e=-r+i 

In analogy to femionic path integrals it is convenient to think of x G M+ as space, j G Z, 
\j\ < t, as time, and t as the number of growth steps. As already anticipated 4> T has 
determinantal moments. This means that for a time-ordered sequence u\ < . . . < u m and 
arbitrary space points x\,... ,x m one has 

m 

(JJ <t>r(xj,Uj)) = det{K {2T) (x t ,u, ; Xj, Uj)}i<i,j< m . (19) 
i=i 

To write down a formula for the kernel K^ 2t ^ we introduce the operators T± on L 2 (R) as 
given by the integral kernels, x, y G M, 

f p-( x -y)/ 2 if t > v 

T +^) = { \lx<\\ 

rp I \ / , if x > y , 

T-(x,y) = l [e _ {v _ x)/2 - lix£y (20) 
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We also introduce the projection P_ onto the half line (—00, 0]. Then, for r > 1, x, y > 0, 

K^\x,2u;y,2v) = {T+T^T {T+TZ 1 ) 7 ■p_(T+ 1 T_) T (T+T_ 

(21) 

in case — t < u < v < t and 

K^\x,2u;y,2v) = (T + T_Y((T + Tl l y P_{T- l T_y - l)(T + T_)~ v (x,y) 

(22) 

in case — t < v < u < t. There are similar formula for odd arguments. 

As a particular case one may consider equal times u, v = 0, i.e. the moments of 
fyirix-, 0). Specializing in (21) results in the equality 

r-l 

(x, 0; y, 0) = £ Lj (z)^ (^e^e"^ 2 , (23) 

x, y > 0, where are the order Laguerre polynomials normalized as 

/ dxe~ x Li(x)Lj(x) = Sij . (24) 
Jo 

Such a kernel is known from the theory of random matrices [23]. There one considers N x 
N complex matrices A where the matrix elements Aij are independent Gaussian random 
variables normalized as (|yljj| 2 ) = 1. Then the eigenvalues A* A have determinantal 
moments with kernel (23). Thus the eigenvalues of A* A have the same joint distribution 
as h e (0,r),£ = 0,... ,-r + l. 

It follows from (19) that the joint distribution function for the top line at times u\ , . . . , u m 
is expressed through the determinant of an operator acting on L 2 (R) (g> C m , i.e. acting on 
m-spinors over M. We choose U\ < . . . <u m and define the integral kernel 

£(x,i;y,j) = if (2T) (x,Ui; y,uj) . (25) 

Let Q be multiplication by the indicator function X[ Sj ,00) i v ), s j > 0. Then 

Prob(/i (ui,r) <si,..., h (u m , r) < s m ) = det(l - g 1 ' 2 ^ 1 ' 2 ) . (26) 

To obtain the asymptotic statistics of h (j, r) thus requires an understanding of the asymp- 
totics of the kernel K^ 2 ^ as r — > 00 at suitable rescaling of the arguments (x, u; y, w). 

4. Extensions 

We discussed the line ensemble for step initial conditions of the TASEP. For large r, the 
dual growth process ho (j, r) has a definite shape, which is bent downwards and has a jump 
of macroscopic size at both endpoints. Thus under edge scaling h (j, r) is governed by 
the stationary Airy process, which has joint distributions of a structure very similar to (26); 
only the kernel K. is replaced by the extended Airy kernel [3,16]. 
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By the same method also general step initial conditions can be handled, to say: we pick 
the left density p_ and the right density p + and require that under the initial conditions for 
the TASEP the rjj's are independent with Prob(?7j = 1) = p_ for j < and Piob(r]j = 
1) = p + for j > 0. In the mapping to the directed polymer the w(i,j)'s have to be 
extended by the row w(j, 0), j > 0, and the column w(0,j), j > 0. The directed polymer 
starts at (0, 0). More precisely the w(i,j) are independent exponentials with the following 
means, 

(w(i, j)) = 1 for i, j > 1 , w(Q, 0) = , 

(w(j,0)) = (l-p+)- 1 , (w(0,j))=pZ 1 for j>0. (27) 

The construction explained in the previous section can be carried through only if 
the w(i,j)'s are independent exponentials with Prob(w(z, j) G [y, y + dy\) = (on + 
exp[— (on + {3j)y]dy, y > 0, oti,(3j > 0. Starting from the left, the weight for 
an up-step of size 8 is then exp[— ctj6] and for a down-step of size S it is exp[— /3jS], 
S > 0. For the step initial conditions, p_ = 1, p + = 0, clearly ~ct = (|, |, . . . ) = [3 . 
In the general case, considered here, the obvious choice is ~b? = (i — p + , |, . . . ) and 
(3 = (p_ — |, |, . . . ), which requires < p+ < ^, | < p_ < 1. Up to an error of order 
1 the TASEP height is then represented by the corresponding height of the dual growth 
process. 

The restrictions on p + , p_ are somewhat surprising, since the directed polymer with 
random waiting times according to (27) is well defined for any choice of p + , p_ . To cover 
the left out domain in densities one has to analytically continue in p + , p_. The details 
tend to be lengthy and we refer to [10] for such a continuation in case of the discrete time 
TASEP. Here we briefly explain stationary initial conditions, p + = p = p_, as the one of 
most physical interest. In this case an + ffo — 0, formally, which means that the weight of 
the very first deposition is 1 independent of its height. Such a weight cannot be normalized 
to 1, but averages of physical quantities can still be computed through a suitable limit 
procedure. As a net result one finds that for joint distribution functions 

Prob p (/io(ui,r) <*!,..., h Q (u m ,T) < s m ) = dct(l - gl /2 ^ s 1/2 ) s =o + 0(1) , 

ds (28) 

compare with (26). Here Prob^ refers to the stationary TASEP with initial Bernoulli mea- 
sure of density p. h (u, r) is the top line of the corresponding line ensemble after r growth 
steps. Q s is the multiplication by the indicator function of the interval [sj + s, oo) and JC 
is the extended Laguerre kernel (21) perturbed by a rank 1 operator which contains the 
information on the density p. Such a structure seems to be novel. 

Step initial conditions lead to boundary terms for the directed polymer. A further natural 
modification is to retain the independent exponentially distributed w(i, j) of mean 1 but to 
alter along the diagonal to Prob(w(i,i) G [y,y + dy}) = r )e~ lv d, y > 0, 7 > 0. This 
results in two well-studied, but unresolved problems: 

(i) For the TASEP the rate for jumps from to 1 is changed from 1 to 7. For 7 <C 1, 
particles pile up behind the slow bond and the average current is reduced. It is conjectured 
that at 7 C = 1 the current reaches its maximal value of j [24]. A more intricate scenario is 
proposed in [25]. 

( ii) On the other hand the directed polymer can be interpreted as an interface in a random 
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potential with a pinning potential along the diagonal {i = j}. For 7 < j c the interface is 
pinned and it unbinds as 7 — > j c . For 7 > j c the directed polymer has n 2 / 3 transverse 
fluctuations [26]. It is conjectured that the unbinding transition occurs at 7 C = 1 [27,28]. 

The exponential distribution of the w(i,j)'s derives from the exponential waiting times 
in the jumps of the TASEP. On the other hand, from the point of view of the directed poly- 
mer any other distribution looks fine, perhaps imposing some moment conditions. Un- 
fortunately the construction of purely entropic line ensembles requires that the w(i, j)'s 
are allowed to be independent with the geometric distribution Prob(w(i, j) = n) = 
(1 — dij)(aij) n , n — 0,1,..., where the <Zjj must be of product form, < ay < 1 
and ciij = bibj. This includes the exponential limit, discussed previously. In the oppo- 
site limit, the Poisson limit, the w(i,j)'s are mostly and once in while take the value 1. 
Thereby one arrives at a Poisson point process in the quadrant {x\ > 0, x^ > 0} with an 
average density, p, of product form, p(x\, X2) = Pit^i)/^^)- The directed polymer lives 
then in the continuum. It moves from Poisson point to Poisson point, linearly interpolat- 
ing, restricted to be increasing in both coordinates and starting from (0, 0) reaching some 
prescribed endpoint {x\ , X2) in the positive quadrant. The corresponding growth process is 
the polynuclear growth (PNG). It has a continuous growth time, t, and the height function 
h(x, t) is over the real line, x 1— > h(x, t) is piecewise constant and takes integer values 
only. The dynamical rules can be read off the directed polymer. Accordingly, an up-step 
moves with speed one to the left and a down-step with speed one to the right. Whenever 
two steps meet they annihilate. In addition, randomly and uniformly in space-time, a spike 
of height one is nucleated on the top of the current height profile. The spike, then immedi- 
ately widens under the deterministic motion. The line ensemble for the PNG is constructed 
as before, namely by using the mass annihilated at line h-i (x, t) as nucleation in the lower 
lying line h-i-i(x, t). The droplet geometry can be imposed by requiring that initially 
h(x, 0) — for x ^ and h(0, 0) = 1 and that nucleation is allowed only on the top of 
the first layer. As for the TASEP, the resulting random field has determinantal moments. 

In general, however, line ensembles do not have to be determinantal. A point in case is 
the PNG half droplet with an extra source of rate 7 at the origin [5]. In the corresponding 
line ensemble the lines are pinned at x = t, i.e. h_i{t,t) — —£, while at x = they 
satisfy a "pairing rule". In particular, the point process {/i_^(0, t),£ = 0, 1, ... } has its 
moments given through a suitable Pfaffian. The top line at x = has in the scaling limit 
for 7 = the same fluctuations as the largest eigenvalue of a GSE random matrix and for 
7 = 1 the one of a GOE random matrix. A further point in case is the PNG model with 
flat initial conditions [6]. While in space-time the process is not determinantal, the point 
process {/i_^(0, t), £ = 0, 1, . . . } has moments given through a Pfaffian. 

Acknowledgements. I benefited a great deal from the collaboration with Patrik Ferrari 
and Michael Prahofer. 
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